Dynamic properties of solitons in the Prenkel-Kontorova Model. Application to 

incommensurate CDW conductors 



Victor M. Burlakov 

Institute of Spectroscopy Russian Academy of Sciences, Troitsk, Moscow 
region, 142092, Russia 



00 
On 



(N 
O 
O 



o 

C/3 



Oh- 



I. INTRODUCTION 



An impact of kink-type solitons on infrared lattice vibrations is studied for incommensurate 
Frenkel-Kontorova model. It is shown that the vibration of particles involved into the kink formation 
is very similar to that in a gap mode around the force constant defect. IR phonon mode intensity 
is found to possess a universal dependence on the system parameters and the kink concentration. 
It is argued that the giant IR peak observed in a number of incommensurate charge density wave 
Q"^ | conductors can be explained in terms of dynamic charge transfer stimulated by kinks. 
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A system of interacting particles in sinusoidal external potential (Frenkel-Kontorova (FK) model JL|) is widely used 
■ for description of a broad variety of physical phenomena, such as statics and dynamics of incommensurate phases 
(see, e.g. ]2|), transport properties in quasi-one dimensional conductors (see Ref. 3 and references therein), adatoms 
diffusion on a metal surface |Q, etc. Peculiar features of the FK model usually explored are related to the kink-like 
solitons. Properties of the kinks have been described in a number of publications |^-|ic[]. The dynamics of the FK 
, model has also been extensively studied, but mostly in relation to dynamics of incommensurate superstructures rather 
7— ( ■ than to the single kink [O-Q. Whereas it is not completely clear yet at what system parameters the single- kink 
effects can be still important. 

Incommensurate charge density wave (CDW) conductors are the physical systems for which the single-kink effects 
can be very important. Vibration and transport properties of these systems attract great interest due to numerous 
peculiarities. The most striking among them are: i) the giant peak of unknown origin in the low frequency infrared (IR) 
spectrum of a number of inorganic CDW conductors such as Kq^MoO^, (TaSe^/ and TaSz [ |l5||l6| ]; ii) nonlinear 
conductivity and noise generation when conducting dc current by these materials. 

The aim of the present study is to investigate an impact of both single kink and kink lattice onto IR active phonon 
spectrum and to specify the range of the model parameters in which its properties can be treated in terms of nearly 
independent kinks rather than in terms of superstructure, associated with the kink lattice. I believe that after 
comparison of the numerical results with the experimental data and proper justification of the model it can serve as 
a bases for understanding of the microscopic nature of the aforementioned peculiar features of CDWs. 

The investigations were performed in two approaches: i) molecular dynamic (MD) simulation was used for the system 
to reach an equilibrium state according to the method proposed in Ref. p7|| , after what all the particles were subjected 
to a small uniform step-like displacement and subsequent vibrations were analyzed via Fourier-transformation; ii) 
eigenvector problem (EVP) was solved in the harmonic approximation to study the vibration spectrum of the system. 
The kinks in this case were taken into account through expansion of the potential energy around particle equilibrium 
positions obtained from MD simulation. 



II. VIBRATION SPECTRUM IN THE PRESENCE OF KINKS 



Let us consider a chain of particles of mass m and charge e with nearest neighbor interaction in the sinusoidal 
external potential V (x) = -jy • cos (2tt ■ -) where a is the potential period. In case of harmonic interparticle 
interaction the motion equation for n-th particle is 

d 2 U n , dU n , „ fnTT TT TT V-a . ( U n \ n ~. 
■ 7 • -TT- + K 2 ■ (2U n - U n -i - U n+ i) + — sm 2tt = e • E(t), (1) 
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where 7 is phenomenological damping and E(t) is external electric field. Let the time dependent position U n of 
the particle can be represented as U n (t) — n ■ a + [7° + S n (t), where is quasistatic variable describing a shift 
of the equilibrium position of the particle with respect to the corresponding potential minimum, S n (t) describes a 
vibration of the particle around the new equilibrium position U®. Then suggesting S n (t) = S n (oj) ■ exp(io; • t) and 
E(t) = Eq ■ exp(itj • t) the Eq. (Q) can be splitted into two equations 

K 2 ■ (2U° U Q n _ x - U Q n+l ) + ^ • sin (2tt • U°) = 0, (2) 



A 2 . . . , , , , e 



S n {to) ■ [V • cos (2tt • U°) - + iuj • 7] + — ■ (2<5„(o;) - <5„_i(o;) - 5 n+1 {io)) = - • Eq. (3) 



Disregarding the trivial case U® = when number of particles N part is equal to the number of potential minima 
Np t, Eq. (0) describes quasistatic kink-like deformation of the chain (due to neglection of the dynamical term we 
restrict our consideration by standing kinks only). Eq. ([}]) describes the particle vibration around the new equilibrium 
position. In the continues limit Eq. (||) reduces to the sine-Gordon equation jL8| with the single-kink solution [ p^[ : 

U®(i) — 2a ■ 7T _1 • arctg |exp ± 2 ^ n ^' a j , Rk = ^ a \J can be considered as the kink radius, i is the kink position. 
Substituting this solution into Eq. (^|) one can obtain the complex susceptibility = -g- • ^<5„(w), where the 

peaks in Im(x(w)) correspond to resonances uj r and Re (5 n (u) r )) corresponds to suitably normalized eigenvector of 
the mode at u) r . 

Fig.l shows the u(k) plot without (a) and with (b-c) kinks in the model system. One can see the smearing of 
the resonances in the presence of single kink (Fig. la) while this smearing is gone in the system with the kink lattice 
(Fig.lc). Instead, the phonon band folding due to the kink lattice and the additional low frequency vibration arise. 
The latter is a so-called phason, which is related to translational motion of kink(s) (domain wall(s)). In case of 
negligible Peierls-Nabarro potential barrier the phason frequency tends to zero. The phason is IR active and looks 
like a steep increase at u> — ► in the optical conductivity spectrum a (u>) = uj ■ Im(x(w)) (see Fig. 2). The high 
frequency peaks in Fig. 2 correspond to phonons, the strongest one being related to in-phase vibration of the particles 
situated near bottom of potential wells. 



III. SIMILARITY BETWEEN KINK AND THE FORCE CONSTANT DEFECT 



The smearing of the vibrational states in the presence of single kink shown in Fig. lb suggests that the kink is acting 
like a point defect. Moreover, the particles involved into the kink formation are almost completely eliminated from the 
high frequency phonon-like normal modes while these particles obviously possess a higher vibration amplitude (local 
vibration density) at low frequencies (see Fig. 3). This is quantitatively illustrated in Fig. 4 where the eigenvectors of 
phason mode and that of IR phonon mode are shown. Accordingly, one could try to describe vibration properties 
of kinks in terms of localized vibrations around force constant defect. That means the original incommensurate 
(Np ar t 7^ Np t) FK model is replaced by commensurate one (or the ordinary harmonic chain of particles with harmonic 
incite and interparticle potentials) and some particle cites possess a defect force constant. The strength of the defect 
has been determined from equation |2(]] 

1 + ~N~ ' £ y + 4X 2 .sin 2 (|) = °' (4) 

what means the zero-frequency gap mode formation in the vicinity of the defect cite. As it is shown in Fig. 4 the 
eigenvector of the gap mode is very close to that of the phason while the eigenvectors of the phonon-like mode nearly 
coincide in both cases. The corresponding spectrum of the ID crystal with the force constant defect is also shown in 
Fig. 2. Note, that the localization length of the gap mode S gap (the halfwidth of the peak shown by dashed line in 
Fig. 4) is equal to Rk/\2 in a wide range of Rk values (see insert in Fig. 4). 

Thus, one may consider the system with kinks as a defect, or impurity crystal taking for the description of its 
vibration properties all the results already known. For instance, it is well understood that S gap is determined basically 
by splitting of the gap mode from the bottom of the optical band lu = W and by the optical bandwidth 2y^2. One 
may argue therefore that the similarity between phason and the gap mode eigenvectors and the phason eigenvector 
itself does not depend on the potential anharmonicity provided that its influence on the above mentioned parameters 
is small enough. It is thought therefore that the obtained results will be applicable for a more realistic interparticle 
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potential too. From the analogy between the kinks and the defects it follows also that the IR phonon mode intensity 
will show a linear decrease versus (rik = is the kink concentration and = \N part — N pot \ is the number of 

kinks) for low kink concentration. This indeed does take place in certain range of Rk values. 

IV. WHEN THE KINK EFFECTS ARE IMPORTANT 

Although the N-kink solution of Eq. (||) is also available J2lj] it is more convenient to approximate it with the 
sum of single-kink solutions. Our MD study of the ground state of a system consisting of 128 particles arranged 
in 128-TVfc potential wells with cyclic boundary conditions showed that even for Nk S> 1 the kink lattice can be 

perfectly described as a sum of the single- kink solutions with Rk ~ 2a^J Namely, for Nk = 8 and K 2 — 4 * V 

(R k heor = 4.0a) the value of R° k xp ~ 3.94a has been obtained. Similar results have been obtained also for the case 

when the number of potential wells exceeds that of the particles. The dipole moment spectrum / (lj) = Im ' J - — 



E 

has been both calculated from (pi) substituting = Yl U®(i) with Rk = R^ 9 and obtained from MD simulation 

»=i 

via fluctuation dissipative approach for various values of r\ = Rk ■ rik- Both approaches agree rather well even at 
very low frequencies although the harmonic approximation obviously fails at oj = 0. Two examples of the particles 
arrangement and corresponding eigenvectors of the IR vibrations are presented in Fig. 5. The eigenvectors for r\ = 0.25 
can be represented as a superposition of the single-kink eigenvectors while for 77 = 0.5 they are looking quite different: 
even those particles which still occupy the potential wells and are not involved into the kinks formation, are strongly 
involved into the characteristic IR vibration (compare (a) and (b) in Fig. 5). It should be pointed out that there is 
no noticeable difference between the commensurate and incommensurate cases (kink lattice period is equal and is 
not equal to an integer number of a, respectively) if the kink concentration is not too high. Otherwise the difference 
manifests itself in a small shift of the zero-frequency peak in Fig. 2 from its position. 

I concentrate here on a question concerned with the intensity of the phonon peaks in Fig. 2 as a function of the 
parameter r\. The EVP approach (Eq. (^) with various values of ^ and rifc) was used for this study. The results are 
presented in Fig. 6. The integrated intensity Is = J I (w) dui of the phonon peaks reveals a universal dependence on 
the parameter r\ S| . It was also found that the eigenvector of the strongest IR vibration obey some sort of scaling 
invariance: the vectors obtained at different rik but for one and the same rj can be transformed to each other by proper 
scaling of a. Note, that the parameter r\ means a volume fraction (in 1-D case) occupied by the kinks and the observed 
decrease in Is at low 77 values can be interpreted as washing out of the high frequency density of states by gap modes 
associated with kinks. At higher 77, when the kinks form real lattice and eventually sinusoidal superstructure due to 
interaction with each other, the decrease in Is oc 77 slows down because the real kink radius can not exceed at least 
one half of the kink lattice period. Indeed, the linear decrease in Is shown in Fig. 6 ends at a cut-off value of 77 ~ 0.4 
(for a = 1) which implies the above mentioned restriction on the kink radius is Rk < 0.4 • k' 1 , (k s = nu) where k s is 
the kink lattice (or superstructure) wave vector measured in units of — . Thus, one can display a range of parameters 

k s ■ J < 0.2 in which the single-kink effects are important, i.e. the system can not be explicitly treated in terms of 

some sinusoidal superstructure related to the kink lattice. Since the IR eigenvectors has been argued to be not very 
sensitive to anharmonicity one might expect this criterion holds for more realistic potentials too. 

V. CDW COLLECTIVE DYNAMICS 

The low frequency excitation spectrum of the CDW ground state have been widely investigated both theoretically 
and experimentally (see, for example, reviews |25j-|25]| and references therein). It is well established that incommen- 
surate CDW ground state is characterized by two specific collective excitations: IR active phase mode and Raman 
active amplitude mode. ]2^| The frequency of the former, u> p , in the incommensurate CDW conductors is of the order 
of 1cm -1 while the amplitude mode frequency uj a is about one-two orders of magnitude higher. These vibrations have 



been observed experimentally in such model CDW conductors as K0.3M0O3, (TaSei)2l and TC1S3 fll6| , |l5| , [27| |32|. 

Besides the phase and the amplitude modes an additional vibration obeying giant IR activity has been observed in 
all the above mentioned compounds ]l5| , p^ |. The frequency of this additional feature in (TaSe4)2l is about 38cm _1 in 
between the phase (~ 1cm -1 Q) and the amplitude (~ 90cm -1 p9| ) mode frequencies. Several explanations have 
been proposed to account for the additional giant IR peak, but microscopic ori gin of this vibration is still not clear 
(see, for instance, the discussions in Refs.l and 15). In phenomenological model j54| the additional giant IR peak was 
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thought to result from a bound collective-mode resonance localized around impurity, but again without emphasis of 
microscopic origin of the model parameters. 

Below it is shown that the giant IR resonance occurs in the incommensurate CDW system even in the absence of 
any impurities provided that the dynamical charge transfer between adjacent CDW periods is taken into account and 
the CDW possesses some kink lattice structure rather than sinusoidal structure. 

Since the lattice deformation coupled to the CDW is much smaller than the crystal lattice constant it is reasonable 
to describe the CDW system within the FK model. The CDW periods in this case are associated with particles of 
mass m and charge e which are placed into sinusoidal external (crystal lattice) potential V (x) — — • cos (2ir ■ ~). 
The interparticle distance in the commensurate phase is accepted to be equal to 2a what means the CDW is formed 
due to dimerization. In case of 2N par t ^ N po t one again obtains incommensurate structure (kink lattice) and the time 
dependent position U n of the particle can be represented as U n (t) — 2 ■ n ■ a + U° + 5 n (t) . 



A. Dynamic charge transfer 



As it has been demonstrated by Itkis and Brill 35 1, spatial redistribution of the charge condensed in CDW takes 
place under action of static electric field. Obviously a characteristic time for the charge redistribution or, in the other 
words, for the charge transfer from one CDW period to another is determined by the amplitude mode frequency 
(~ 90cm- 1 Hi). Therefore, in the case of the giant peak frequency the adiabatic condition is fulfilled. 

To take into account the charge transfer contribution to the IR intensity of any mode let us suppose that the particle 
charge in our model is determined as 



= e • (1 +0 ■ (S n+ i(t) - <Wi))), 



(5) 



what means that the charge is transferred from the region of local compression of the CDW to a region of local 
dilatation. The factor /3 determines the fraction of the particle charge transferred during vibration. The dipole 
moment Pit) is determined as a sum of the part related to the particles displacement Po(t) and the part, related to 
the charge transfer between adjacent unit cells P c t(t) 



P(t) = P (t) + P ct (t) = 
e ■ [S n (t) + 0-(U°- • (S n+1 (t) - Sn-i® - S n {t) + *„_ a (t))] 



(6) 



In commensurate phase U® — U®_ 1 = a for all n and the charge transfer dipole moment P c t{t) vanishes according to 
(||). In incommensurate phase the P c t{t) value can be rather high. It will be shown that the charge transfer effect is 
essentially determined by kink- related disturbance of the periodicity in the particle arrangement. Using the criterion 
of importance of single-kink effects obtained in the preceding section one can examine if the kinks are important 
for description of the charge density wave conductor (TaSei)2l ■ The superstructure wave vector in this system is 
k s ~ 0.085 p6|], \JVjm can be associated with the giant IR peak frequency lu ~ 0.005 eV 37 and y/K-z can be 



estimated from the phason dispersion y/ Ki j m < 0.001 eV [ pSfl . Thus, one obtains k s ■ J <c 1 what implies that 
the kink effects can be important in this compound. 

Fig. 7a shows the fragment of MD simulation of arrangement of 128 particles over 264 potential wells, i.e. the 
CDW with superstructure. The 51-th particle (shown by arrow in the figure) is pinned. The conductivity spectra 



I (uj) — uj ■ Im 



obtained from MD simulation are shown in Fig. 8 for both pinned and depinned system. The 



features of the interest are the phase mode (PM) and the peak of CT mode (charge transfer mode), marked by PM- 
and CT-arrows in Fig. 8, respectively. The latter peak is genetically related to the vibration with the wave vector 
equal to that of the superstructure. The corresponding eigenvectors are shown in Fig. 7b. Taking into account that 
the CDW internal deformation can be adiabatically accompanied by charge redistribution one obtains that the CT 
peak acquires the giant IR intensity. The conductivity spectra in which the charge transfer effect has been taken into 
account according to Eqs (|5|) and (||) are shown in Fig. 8 by symbols (depinned chain) and thin solid line (pinned 
chain). The phase mode intensity is almost independent on the charge transfer effect while the CT mode intensity 
increases several orders in magnitude. 



Fig. 9 shows that the CT mode intensity Im 



Eg 



with charge transfer contribution decreases with the 



increase of the parameter 4au/ and possesses a universal dependence regardless the superstructure wave vector. 
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Physically this feature can be understood in the following way. The charge transfer dipole moment consists of a sum 
of elementary dipole moments resulting from the charge transfer over the inter-kink distance. The longer is the latter 
the higher is the elementary dipole moment, but the smaller is the number of these dipole moments. Thus, the total 

charge transfer dipole moment does not depend on the kink concentration (probably unless 4a^/ ^ < a/rife). On the 
other hand, this dipole moment strongly depends on the kink-mcdiatcd distortion of the chain. The latter decreases 
with the increase of the kink radius resulting in the observed decrease in CT mode intensity in Fig. 9. 

It was experimentally proved that the giant IR peak intensity in {TaSe^I increases with the increase of the sample 
temperature p9| ] . This unusual feature can be naturally explained in terms of the charge transfer effect discussed 
above. Indeed, as it is clear from Fig. 9 the integrated intensity of the CT mode as a function of the system parameters 
can be approximated as (see dashed line in Fig. 9) 



C -e 



2 



n -2.25 

K 2 



(7) 



where Co is constant. The interparticle force constant, obviously, depends on the particle charge K% — e 2 K 2 since 
the particles are associated with the charges condensed in CDW. The (0) can be then rewritten as 



I ~ C n • e 



-0.25 




-2.25 



(8) 



Thus, despite the decrease of the particle charge (the CDW amplitude) the CT mode integrated intensity increases 
upon approaching T p ~ 261K from below! Due to short range order the particle charge (CDW amplitude) remains 
finite even at very high temperatures and accounts for the high CT mode intensity well above T p . 

Note one more peculiarity of the CT peak. It does not depend on the number of particles in the coherent CDW 
domain or, in the other words, on the effective mass of CDW condensate. The latter can explain why the corresponding 
frequency has nearly the same value in such different compounds as K0.3M0O3 and (TaSe^-f ]3^ , |33[ |. 

VI. SUMMARY 

The kink-like solitons in the incommensurate Frenkel-Kontorova model are investigated regarding to their impact on 
the vibration spectrum. It is found that the IR phonon intensity possesses universal dependence on the product of the 
kink radius and the kink concentration suggesting some sort of scaling invariance for the corresponding eigenvector. 
The model accounting for the giant IR peak in the incommensurate inorganic CDW conductors is proposed. It is shown 
the giant IR peak is related to the fundamental vibration with the wave vector equal to that of the superstructure 
and the giant IR intensity is caused by dynamical charge transfer accompanying the CDW internal motion. 
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FIGURE CAPTURES 



Fig.l. Dispersion of vibration band obtained via EVP solution (system of linear equations (|3|)) for the FK model 
containing 64 particles arranged over: (a) 64 potential wells (no kinks); (b) 63 potential wells (one kink); (c) 56 
potential wells (four kinks). The dark regions correspond to the higher vibration amplitude of particles excited by 
external field E = E cos(fc • n) ■ cos(wi). 

Fig. 2. Conductivity spectra of the FK model with one kink. (1) is the calculated spectrum; (3) is that obtained by 
MD simulation (32 particles with cyclic boundary conditions) for K2 = 4 • V, yV = 72 arb.un; (2) is the spectrum 
corresponding to the force constant defect AU = —4.1231 • V (see text), in the position of the particle no. 16. 

Fig. 3. MD study of local density of states distribution over the particles in the commensurate (a) and containing 
one kink (b) FK model. The particles were initially subjected to random displacements and the temporal evolution 
of the spatial harmonics was analysed via Fourier-transformation. The dark regions correspond to the maxima in the 
Fourier spectrum. 

Fig. 4. The kink and the phonon (the strongest peak in Fig.l) eigenvectors obtained as described in the caption to 



Fig. 5. (a) The particle arrangement in the FK model of 128 particles (shown by symbols) in 120 potential wells 
(solid line), (b) The eigenvectors of the kink-(l,3) and phonon-like (2,4) modes. K2 = 4 • V (solid symbols) and 
K 2 = 16 • V (open symbols), \JV = 72 arb.un. 

Fig. 6. Integrated intensity of the phonon-like modes upon r\ = ■ n^ calculated using (0) for FK model of 128 
particles arranged over: (1) 112 potential wells (n^ = 1/8); (2) 120 potential wells (rife = 1/16); and (3) 124 potential 
wells (n k — 1/32). 

Fig. 7. (a) Fragment of particle arrangement obtained via MD simulation in the FK model containing 128 particles 
arranged over 272 potential wells (nk = 1/16, see text) for K2 = 16 • V. The particle No. 51 is pinned by an extra 
local potential, (b) Eigenvectors of phason mode without pinning (dotted line) and with pinning (thick solid line), 
and those of CT-modc (sec Fig. 2) shown by thin solid line for both pinned and depinned chain. 



Fig. 8. IR conductivity spectra of the FK model (see caption to Fig. 7), calculated using Eqs. (3)-(5) for 4y -y- = 8 

(a = 1). Thin solid line is for (3 — and thick solid line is for (3 = 30 in case of pinned chain. Dashed line is for 
(3 = and symbols are for f3 — 30 in case of depinned chain. PM are the phase modes and CT are the modes which 
intensity may contain the charge transfer contribution. For (3 = 30 the 0.03 • e of the particle charge is transferred 
during the CT mode vibration while for (3 = it is 0. 

Fig. 9. Dependence of the CT mode integrated intensity on the potential parameters of the FK model containing 128 
particles arranged over 264 minima (rife = 1/32): (1) is for f3 = 30 and (4) is for [3 = 0; over 272 minima (nk = 1/16): 
(2) is for (3 — 30 and (5) is for (3 — 0. (3) is the dependence y = 0.45/x 2 25 . Arrow shows the parameters for the 
spectra presented in Fig. 8. 



Fig.l. By symbols in the insert is shown the dependence of the gap mode radius upon kink radius Rk = 2 
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